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We study the Luttinger-liquid parameter Kp of the Hubbard chain and the Hubbard lad- 
der models by the ordinary perturbation method combined with the Luttinger-liquid relation. 
According to the Luttinger-liquid relation, the critical exponent Kp is related to the charge 
susceptibility Xc and the Drude weight D by Kp = ^{nXcD)^^^ . By calculating these quanti- 
ties with the perturbation method, we obtain Kp at the first-order analytically and up to the 
second-order numerically. We compare these results with results of the Bethe ansatz for the 
Hubbard chain and that of the numerical diagonalization for the Hubbard ladder. It shows that 
the perturbation calculation of Kp is reliable in the weak coupling regime. 
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§1. Introduction 

Low-dimensional strongly correlated electron systems 
have been investigated extensively due to the possible 
relevance to high- Tc superconductivity. In particular, 
One-dimensional(lD) Hubbard-like models provide an 
important ground for understanding electron-correlation 
effects. It is well known that a single Hubbard chain 
has generic properties of Luttinger liquid state with gap- 
less spin and charge modes for repulsive interaction. For 
attractive interaction, a superconducting state charac- 
terized by idpHunant paring correlations with a spin gap 
is reahzed.tl'BEl' 

Recently, weak coupling theory has been applied to 
the problem of ladder models, which are interesting as 
a first step towards 2D systems and may be relevant for 
some materials. The theory reveals that the systems 
remain typical jiqiil- Fermi-liquid properties as ID elec- 
tronic systemsMcl'Q^ At half-filling, they are Mott in- 
sulators, exhibiting gaps to all excitations. Upon dop- 
ing, the gaps survive except one gapless charge mode 
and a superconducting (SC) paring correlation charac- 
terized slowly power-low decay appears. According to 
the Luttinger-liquid theory, critical exponents of vari- 
ous types of correlatiDn. lunctions are determined by a 
single parameter iirp.El'EH3) It is predicted that the SC 
correlation function of ladder systems is dominant for 
Kp > 0.5 at hghtly doping. It decays as ~ r ^^^p , 
whereas the "Akp" charge density wave (CDW) correla- 
tion function decays as ^ r~^^p . The spin density wave 
(SDW) and "|2^" CDW correlation functions decay ex- 
ponentially.oB' EP 

In spite of the good understanding of the single Hub- 
bard chain which is exactly solved by the Bethe ansatz 
method, Hubbard ladder models are much less known. 
In fact, systematic treatment of Kp as a function of in- 
teraction U is not yet obtained. In this work we propose 
a simple method for calculating the Luttinger-liquid pa- 
rameter Kp. It is the first and second-order perturba- 
tion expansion with respect to U combined with the 



Luttinger-liquid relation. These results are compared 
with the Bethe ansatz and the numerical diagonalization 
results. 

§2. Perturbation expansion of Kp 

In the Luttinger liquid theory, some relations 
have been established _as universal relations in one- 
dimensional models E'D'cP The critical exponent Kp is re- 
lated to the charge susceptibility x and the Drude weight 
Dhy 
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where N is the number of lattice site, n is electron den- 
sity and Eg is the total energy of the.^ ground state as 
a function of n and a magnetic flux (/)Jj' If we calculate 
the ground state energy in some way, we can obtain Kp 
through X and D. 

We apply the ordinary perturbation expansion in pow- 
ers of U to E„ 



Eg = Eo + ei + e2 ■ 



where Eq is the ground state energy of the non- 
interacting system and Ci (?=1,2,....) is the i-ih order 
correction of Eg. According to the perturbation expan- 
sion of Eg, x^^ and D are given by 



X ^ ^Xa^+xiU+X2U\..., D = Do+diU- 

where Xo^ and Do{— ^/tt) are the inverse charge sus- 
ceptibility and the Drude weight of the non-interacting 
system respectively. Coefficients Xi and di (i=l,2,....) 
are the i-th order corrections of x~^ and Z?, which 
are determined by the second differential coefficients of 
Ci with respect to n and (j) respectively. In the Hub- 
bard model, the first-order term ei is easily obtained 
as ei = UJ2i < "iT X "■ii UNn'^/A. It leads 
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xi = \/2 and di — 0, where ei is independent of the flux first-order expansion of U 



Substituting these values into eq. (1), we obtain Kp 
within the first perturbation expansion: 



(5) 



(3) 



It shows that Kp generaUy decreases with increasing the 
repulsive interaction U in the weak coupling limit. It 
is noted that this approximation is equivalent to the 
Hartree-Fock (HF) approximation. The second-order 
term of the ground state energy 62 is also easily obtained 
as 

g V- /(fci)/(fc2)(l - /(fci + g))(l - /(fc2 - q)) 

N ^ eih) + sik2) - e{h + q) + sik2 - q) 



where e(fc) is the non-interacting band and /(fc) — 
6{kF — \k\). Using 62, we can calculate Kp up to order 
C/2 by 



l/Klc^l + ^U+{x2-^d2)xoU^ 



(4) 



with 



X2 



1 9^62 



d2 = 



In the following section, we will estimate the above inte- 
gral and obtain X2 and d2 numerically. 




Fig. 1. The coefficients of the second-order terms X2 and 6,2 as a 
function of n. The broken line represents {x2 — rf2)- The inset 
shows the ground state energy of the system as a function of U 
by the first-order and the second-order perturbation calculation 
with the exact result at quarter-filling. 



§3. 



Hubbard chain and Hubbard ladder 



At first, we examine the Hubbard chain model whose 
non-interacting band is given by e{k) — —2tcosk. In 
this case, we have Xo^ ~ T^tsinkp and Kp up to the 



where kp = is the Fermi wave number. This expres- 
sion of Kp is exactly equal to the result of bosonization 
method for the Hubbard chain. E|) 

To estimate 62, we use N = 200 and 400 sites systems. 
We confirm that the size dependence of X2 and c?2 are 
very small and negligible. In Fig.l, we show X2 and d2 
as a function of n. These values seem to diverge at the 
limit ri — > 1. It might reflect the insulator transition 
of the Hubbard chain at half-filling. We also show the 
ground state energy Eg /N as a. function of U at quarter- 
filling. The result of Eg/N is consistent with the exact 
result in the weak coupling region. In Figs. 2(a) and (b). 





Fig. 2. Kp as a function of n by (a) the first-order and (b) the 
second-order perturbation expansion (b) with the exact result 
of the Bethe ansatz. The solid lines represent the results of the 
Bethe ansatz and the broken lines are that of the perturbation 
expansion. 



we show Kp by the perturbation expansion with the ex- 
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act result of the Bethe ansatz.ll3'tll' They show that the 
result of the perturbation expansion is consistent with 
the exact result in the weak coupling region. It also 
indicates that the second-order perturbation calculation 
agrees with the exact solution more than the first-order 
calculation in the weak coupling regime. In the strong 
coupling regime, Kp of the first-order calculation seems 
to be close that of the exact solution at n ^ 0.4. How- 
ever, it may be an accidental. 

Next, we consider the Hubbard ladder Hamiltonian 

(6) 

where c\ ^ ^ stands for a creation operator of an electron 
with spin a at site (i, a) and U is the on-site interaction. 
Here, a{= 1,2) denotes legs and i is rung. In this case, 
non-interacting band is written as 

£^(fc) = —2ti cos k ± 

where £+(fc)(e~(fc)) represents the upper (lower) band 
and k is the wave vector. If we define kp_{kp^) as 
a Fermi point in the lower (upper) band, we find that 
e^{kF-) ~ e^{kF+) with kp_ + kF+ = kp = mr. The 
differential coefficient of e^{kp^) with respect to n leads 
the inverse charge susceptibility Xo ^- After a bit of cal- 
culation, we get 

_i ^ . , ,coskp{l-coiikp) + {tr/tiY -'^'"'^'^ 
Xo =-t,7rsmfcF{ 



On the other hand, for n < 0.5, the values X2 — f ^2 
are negative. In particular, the absolute values are very 
large near n = 0.5. It indicates that Kp is enhanced 
by the second-order term. In contrast to the case for 
n > 0.5, the existence of the upper band seems to pro- 
duce an effective attraction in the second-order correc- 
tion. This result may consist with an enhancement of 
the SC correla-tion shown by the numerical diagonaliza- 
tion methodO Figure 4 indicates Kp as a function of n 
at U/tr = 1.0 by the first- and the second-order pertur- 
bation calculation. 
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Fig. 3. The coefficients X2 and d2 as a function of n for the Hub- 
bard ladder modeh The broken hnes represent {x2 — ^d2). 



(1 — coskpY 



X {{tr/UY 



, [tr/tl 



sin kp 



1 — cos kp 



-1/2 



Substituting Xo to eq.(3), we have Kp of the Hubbard 
ladder model axialytically within the first-order pertur- 
bation method.Ej' 

When n is smaller than Uc which is determined by 
£+(0) — e^(^^^), electrons are filled only in the lower 
band. The density of state of the lower band is a 
half of that of the chain model. Then, xo is given by 
l/(27rsinfci?) which is a half of that of the chain model. 
Within the first-order perturbation calculation, the up- 
per band is irrelevant and the correction of Kp becomes 
a half of that of the chain model. 

As well as the Hubbard chain model, we estimate the 
second-order terms 62, X2 and c?2 numerically.li3 In Fig. 3, 
we show X2, d2 and X2 — -1^2 as a function of n at 
tr/tl = 1.0. For n > ric = 0.5, the values of X2 seem 
to be small, but the absolute values of c?2 are large. It 
shows that the values X2 — jd2 are positive and large. 
It indicates that the second-order term of the repulsion 
U reduces Kp as well as the first-order term. This re- 
sult does not contradict the behavior of the SC corre- 
lation function obtained by density matrix renopnaJiza- 
tion group and quantum Monte Carlo methods O'Ej'Ej^ 
Although system sizes used in these numerical works are 
too small to determine Kp precisely, we find that a rough 
estimate indicates 0.5 ^ Kp < 1.0. 




Fig. 4. Kp as a function of n for the Hubbard ladder model by 
the first-order (the broken lines) and the second-order (the solid 
lines) perturbation expansion. 



To clarify the validity of the perturbation expansion of 
Kp, we examine a finite size system of the ladder model. 
We numerically diagonalize the Hamiltonian of 14 sites 
(7 unit cells) system by using the Lanczos algorithm. 
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We use the periodic boundary condition for A^e = 8, the 
Moebius boundary condition for A^e = 10 and the an- 
tiperiodic boundary condiiion for A^e = 12, where is 
the total electron number .tS This choice of the boundary 
condition gives either fuUy occupied or empty single par- 
ticle orbitals ( closed shell ) and removes accidental de- 
generacy in the non-interacting case. The uniform charge 
susceptibility Xc and the Drude weight D is calculaigd 
from the ground state energy with the usual method.Lr 



0.6. 




U 



Fig. 5. Kp of the finite size system for the Hubbard ladder as a 
function of U. The soUd circles represent the result of the nu- 
merical diagonalization for 7 rungs system at tr/ti = 0.6. The 
broken line represents the result of the first-order perturbation 
expansion and the dashed line is that of the second-order per- 
turbation expansion. The inset shows the ground state energy of 
the system as a function of U by the first-order and the second- 
order perturbation calculations with the result of the numerical 
diagonalization. 



In Fig. 5, we show Kp and the ground state energy 
Eg/N of the finite size system with the result of the nu- 
merical diagonalization method. In the weak coupling 
regime, the results of the perturbation approximation 
and the numerical diagonalization are in agreement with 
each other. It also shows that the second-order pertur- 
bation is better than the first-order in the weak coupling 
regime. 

§4. Summary and discussion 

In this work we examine the Luttinger-liquid param- 
eter Kp of the Hubbard chain and the Hubbard lad- 
der models by the ordinary perturbation method com- 
bined with the Luttinger-liquid relation. According to 
the Luttinger-liquid relation, we obtain Kp at the first- 
order analytically and up to the second-order numeri- 
cally. Comparing Kp with the exact result of Bethe 
ansatz and that of the numerical diagonalization method, 
we show that the analysis of perturbation method is re- 
liable in the weak coupling region. 

Generally speaking, the validity of the perturbation 



lExpansion is not always obvious in ID electron systems. 
\}3 However, it has been analytically shown that the ordi- 
nary perturbation expansion of the ground state energy 
of the Hubbard chain agrees with the expanspp-of the 
exact Bethe ansatz solution at half-filling.HEl'E3 Al- 
though the convergence radius is zero, the perturbation 
is possible as an asymptotic expansion. 

Away from half-filling, the integral equations of the 
Bethe ansatz solution are not analytically solved. How- 
ever, we confirm that 62 consists with the second-order 
term of U in the Bethe ansatz solution numerically. It 
suggests that the perturbation expansion inxipwers of U 
is also possible as an asymptotic expansion.c3 
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